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Abstract 

In their comment arXiv:1212.3831vl Y. Zhou and M.W. Wu claim that the fundamental trans- 
port equation relating the current to the transmission function, used by us and in fact by numerous 
other researchers, is invalid for extended systems and should be corrected. They provide a "cor- 
rect" new formula for transport in extended systems. This would be indeed a surprising new aspect 
of quantum transport theory. Here we show mathematically, however, that the new formula is a 
misconception resulting from adding an energy and momentum dependent function that has to 
vanish due to fundamental reasons. Results and conclusions stemming from adding this function 
are irrelevant. The known established formulas for quantum transport are consistent with each 
others under the well-documented conditions. 
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In their comment [T] arXiv:1212.3831vl Y. Zhou and M.W. Wu argue that the steady- 
state current / in an extended quantum system connected to two leads L and R evaluated 
according to [2H7] 

I = Ilb = -^ j detr {TlG+TrG-) (1) 
is not consistent with / when calculated with the other established formula [6H8] 

/oo 
de tr [ihTL - /ijr^)(G'+ - G') + (Tl - Tj,)G<] (2) 
■oo 
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(for finite systems they claim I^b = Imw)- We use the standard notations where fa = 
fa{t — fia)', <y = L, (R) is the distribution function on the left (right) lead with the chemical 
potential /i^ (fiR) and Ta = i (S+ — where is the respective selfenergy. The super- 
script + (— ) stands for retarded (advanced) quantities. The operator G"^ satisfy a Dyson 
equation that can be written as 
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- (s! + s^) (3) 



hereby gdz) = {z — H^^^ is the resolvent of He-, which is a hermitian Hamiltonian de- 
scribing the central region and gf{z) = gc{e ± irj) where is a small positive real number 
taken to zero after performing the trace and the energy integration in eq.(T][2). G^ = 
iG+ihTL + fRTR) G-. 

Zhou and Wu claim that our previous results [9], [10] calculated with eq.([T]), and for that 
matter the results of all other researchers employing the same approach for an extended 
system, lack scientific ground. Let us show mathematically that the inconsistencies found 
by Zhou and Wu when using eq.([T]) vs. eq.([2]) are self-made and indeed eq.([T]) and eq.([2]) 
should yield consistent results independent of whether the system is finite or extended (in 
the sense introduced by Zhou and Wu in Ref.[T]). 

For clarity let us work in a representation free manner and write the operator equation 
{G+ - G-) = G+ (g—^ - G+^^j G- which readily yields 

G+-G- = G+ [{e -iri- H,) - {e + i7] - H,)] G" - iG+{rL + Tn)G- (4) 
= -2t7]G+G- - iG+iTL + rR)G-. (5) 

The trace of this equation is to be compared with the eq.(5) in Ref.p. Zhou and Wu 
claimed that the first term in Eq. ^ vanishes ... in the finite system, in consistence with the 
previous literature. However, in the infinite system discussed by Ding et al. the situation 



becomes totally different. Hence they derive a new correct formula for extended systems by 
taking the first term of Eq. ([s]) into account and construct a way to make it finite. 

As a matter of fact for a finite Fq,, mathematically the trace of the first term —2iriG^G~ 
of Eq. ([5| has to vanish always and in any basis when 77 — )■ 0, for in this case has neither 
isolated poles nor a branch cut for 77 — )■ 0, i.e. when approaching the real energy axis. This 
is evident from the structure of G^ (cf. [s]). For Fq, — )■ 0, the trace over —2irjG~^G~ yields 
for ?7 — )■ indeed the spectral density of the system, however the current Imw vanishes as 
clear from eq.([2]) (because F^ = 0). This means in turn that introducing somehow a finite 
trace over —2iriG^G^ for r/ — )■ even though Fq — is finite, amounts to a change of the 
spectral density and raises thus the question of the charge conservation in this case. 

We infer mathematically thus that for steady-state transport, i.e. when F^ is finite, the 
first term of Eq. ^ 

-2i7]G+G- = G+ [(e -i7]- H,) - (e + ir/ - H^)] G' 

plays no role. Nonetheless, Zhou and Wu argue that the term —2iriG^G~ should be finite 
because in eq.Q one may write 



[(e - ^r/ - H,) - {e + zr] - H,)] G' = G+ [g; ' - gt '] G". (6) 

While one may do such a manipulation the contribution of this term to the current remains 
of course zero. Zhou and Wu generate nonetheless a finite contribution of this term in the 
limit of —7- by assuming erroneously that even if Fq, is finite (to allow for transport) 
G^ and g"^ are diagonal in a common basis, which is a contradiction as explicitly shown 
above. A common basis allows namely to insert for gf in eq.([6]) the respective spectral 
representation in the absence of the lead when taking the trace in this basis over {G^ — G^) 
(cf. [11] below). From this simple analysis it is obvious that the physical meaning of eqs.(5- 
8) in Ref.[T] is unspecified. Clearly, replacing arbitrarily a vanishing term by an energy and 
momentum dependent function may lead to a series of conclusions that are at variance with 
known results, including the statement that eq. ([T]) is not applicable for an extended system. 

The matter of fact however, for a finite Fq only the second term of eq.(|5]) contributes 
when taking the trace in eqs. ( [Ip| ) and for — > 0. It is straightforward to show that by 



inserting eq.([5]) into eq.([2]) one retrieves the established result given by eq.([T| that we and 



others use for the calculation of the steady-state current, i.e. 

/ = r / detr{ih- fn) [TLG+TnG- + r^G+rz.G'-] } 



(7) 



Hence, as far as the system size is concerned, as introduced by Zhou and Wu and we 
only focus on this issue here, one may use eq.Q or eq.Q and finds Ilb = Imw- More 
importantly any effects on the transport based on a finite first term in eq. § should be 
considered artificial and resulting from some uncontrolled approximations. Based on their 
"correct" formula Zhou and Wu go even a step further in their conclusions and state that 
works for an extended system employing eq. ([T]) are incorrect and lack scientific ground. This 
statement is clearly a consequence of a self-made finite term that for a finite current should 
be in fact zero due to fundamental reasons. 
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FIG. 1. (color online) TMR in a graphene monolayer in contact with two metallic ferromagnetic 
leads as a function of the applied bias V, calculated according to eq.Q or eq. ^ while correctly 
neglecting the first term of eq.Q (as done in our previous work [9|). The dots are the results 
derived by setting the energy width D as the unit of energy. The solid line is the result by using 
the broadening function Tq as the unit of energy. D/Tq = 12 and the spin polarization of the two 
metallic ferromagnetic leads is 40%. The leads are assumed to be of the same material. 
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FIG. 2. (color online) as in Fig. [T| however the TMR is calculated by using new suggestion by 
Zhou and Wu in Ref.^. The parameters are the same as those in Fig. [l| 



Based on their new, allegedly "correct" formula Zhou and Wu raise some issues concerning 
our results [HI HO] , in particular those for the Tunnelmagnetoresistance (TMR) of a graphene 
monolayer contacted to metallic ferromagnetic leads that we obtained on the basis of eq.([T]). 
Since our results do not agree with their calculations based on their own constructed formula 
they claim that the reason is due to a "wrong" energy cutoff D that we use in our calculation. 
It should be noted, that as shown by Zhou and Wu the spurious term they include in their 
formula contains D in a non-trivial manner. In their comment p] they make several claims 
based on the dependence of their calculated current on D. While this point is somehow 
technical, a clarification might be useful to avoid a misunderstanding of the meaning of D. 
We choose D as to ensure the conservation of the number of states in the Brillouin zone 
(after linearizing the spectrum, cf. Ref.[T2]). Calculating D accordingly one arrives at the 
Green function given in our works. 

In fact as shown in Figs,T][2, we can perform the calculations using D or Fq (Fq, are 
assumed to be momentum and energy independent) as the energy scale and arrive at the 



same behaviour of the TMR. Proceeding as suggested in the comment arXiv: 1212. 383 Ivl by 



Y. Zhou and M.W. Wu one arrives at opposite physical behaviour of TMR, i.e. a zero TMR 
at zero bias instead of a peak. Indeed, it is straightforward to show analytically, that this 
behaviour is a direct consequence of assuming, as done by Zhou and Wu, a finite first term 
in eq.Q that derives from the spectral density of graphene. In addition, in their comment 
Zhou and Wu show the result of their "correct" formula for the conductance which exhibits 
basically a linear increase with increasing the bias. This spurious behaviour can be shown 
mathematically to be the direct consequence of including erroneously a finite first term in 
eq.([5| which for small bias dominates near the Dirac point, i.e. the wrong term provides the 
major contribution. 

In Summary, eq.([T]) and eq.([2]) are valid irrespective of whether the system is finite or 
extended. The claims of Y. Zhou and M.W. Wu in larXiv:1212.3831vT! are the result of a 
fabricated finite energy and momentum-dependent term that should vanish due to funda- 
mental reasons if the current is finite. Established approaches to quantum transport are 
consistent within the well-documented limits. 
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